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Abstract. We have carried out an analysis of singularities in Kohn variational 
calculations for low energy (e+ — H2) elastic scattering. Provided that a sufficiently 
accurate trial wavefunction is used, we argue that our implementation of the Kohn 
variational principle necessarily gives rise to singularities which are not spurious. We 
propose two approaches for optimizing a free parameter of the trial wavefunction in 
order to avoid anomalous behaviour in scattering phase shift calculations, the first 
of which is based on the existence of such singularities. The second approach is a 
more conventional optimization of the generalized Kohn method. Close agreement is 
observed between the results of the two optimization schemes; further, they give results 
which arc seen to be effectively equivalent to those obtained with the complex Kohn 
method. The advantage of the first optimization scheme is that it docs not require an 
explicit solution of the Kohn equations to be found. We give examples of anomalies 
which cannot be avoided using either optimization scheme but show that it is possible 
to avoid these anomalies by considering variations in the nonlinear parameters of the 
trial function. 



PACS numbers: 02.10.Yn, 34.80.Uv 



Submitted to: J. Phys. A: Math. Gen. 



Anomaly-free singularities in the Kohn method 



2 



1. Introduction 

Despite the absence of an explicit minimization principle, variational methods have been 
used successfully in many problems of quantum scattering theory. Such calculations 
typically exploit a stationary principle in order to obtain an accurate description of 
scattering processes. The Kohn variational method [1] has been applied extensively to 
problems in electron-atom [2] and electron-molecule [3l H] scattering, as well as to the 
scattering of positrons, e+, by atoms [El El [7] and molecules [8], [HI ^0\. It has been 
widely documented, however, that matrix equations derived from the Kohn variational 
principle are inherently susceptible to spurious singularities. These singularities were 
discussed first by Schwartz [HI [12] and have subsequently attracted considerable 
attention [131 [TH [151 [16] . the region of these singularities, results of Kohn calculations 
can be anomalous. 

Although sharing characteristics similar to those exhibited by scattering resonances 
[T5] , Schwartz singularities are nonphysical and arise only because the trial wavefunction, 
used in Kohn calculations to represent scattering, is inexact [13]. For projectiles of a 
given incident energy, anomalous results are confined to particular formulations of the 
trial wavefunction and can, in principle, be mitigated by a small change in boundary 
conditions or some other parameter. It has also been shown [H [T71 [18] that the use of 
a complex-valued trial wavefunction avoids anomalous behaviour except in exceptional 
circumstances. Alternative versions of the Kohn method have been developed in terms 
of a Feshbach projection operator formalism [19] and have been found [20] to give 
anomaly-free results. 

In this article we will discuss our investigations of Schwartz-type anomalies for 
generalized Kohn calculations involving the elastic scattering of positrons by molecular 
hydrogen, H2. We will find that our choice of trial wavefunction contains a free 
parameter that can be varied in such a way as to produce singularities which are 
legitimate in the context of the scattering theory and which do not give rise to anomalous 
results. Indeed, these singularities can be used to formulate an optimization scheme 
for choosing the free parameter so as to automatically avoid anomalous behaviour 
in calculations of the scattering phase shift. The novelty of determining the phase 
shift in this way is that an explicit solution of the linear system of Kohn equations is 
not required. We will also develop an alternative optimization and show that the two 
schemes give results in close agreement. Further, the results obtained will be seen to be 
in excellent agreement at all positron energies with those determined via the complex 
Kohn method. 

We will give examples of anomalous behaviour which cannot be avoided with either 
optimization, and show that the same anomalies appear in our application of the complex 
Kohn method. We will discuss circumstances under which these anomalies might occur. 
We will show also that such results are nonphysical by considering small changes in the 
nonlinear parameters of the trial wavefunction. 

Our investigations of singular behaviour have been carried out as part of a wider 
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study on (e+ — H2) scattering and annihilation using extremely flexible wavef unctions. 
Our ability to recognize clearly and analyze the anomalous behaviour is as good for 
this system as it would be for a simpler model system, with the advantage that our 
calculations can be used to provide meaningful and physically relevant results [TO] . 

2. Theory 

2.1. The generalized Kohn variational method 

The Kohn variational method is used to calculate approximations to exact scattering 
wavef unctions. Determining an approximation, \i/t, allows a variational estimate, rj^, of 
the scattering phase shift to be calculated, the error in which is of second order in the 
error of from the exact scattering wavefunction, [21]. The standard approach in 
Kohn calculations is to assume an overall form for \i/t that depends linearly on a set of 
unknown parameters, optimal values for which are then determined by the application 
of a stationary principle. 

In our investigations of anomalous behaviour in Kohn calculations for (e"^ — H2) 
scattering, we have studied the lowest partial wave of S+ symmetry. This partial 
wave has been shown [8] to be the only signiflcant contributor to scattering processes 
for incident positron energies below 2 eV. The flrst signiflcant inelastic channel is 
positronium formation which has a threshold at 8.63 eV [21]. Although we will here 
consider positron energies higher than these thresholds, it is not our intention to provide 
a comprehensive physical treatment of the scattering problem taking higher partial 
waves and inelastic processes into account. The purpose of the present study is to give 
a correct and, as far as possible, anomaly-free treatment of the lowest partial wave. It is 
important to examine the single channel case as accurately as possible as a preliminary 
for more sophisticated calculations. By not taking into account additional channels, it is 
possible that anomalous behaviour could occur due to physical inaccuracies in the trial 
wavefunction at higher energies. However, we will demonstrate that all of the anomalies 
in our results ultimately can be attributed to purely numerical effects. 

We have initially used a trial wavefunction having the same general form as 
described in our earlier calculations [9], 

M 

^t = {S + atC + poXo)i^G + ^PiXu (1) 

i=l 

where 
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for some phase parameter, r G [0, tt), with 
N 

S=j^^ri[c{\s-1)], (3) 
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and 

C = cos [c (A3 - 1)] {1 - exp [-7 (A3 - 1)]}. (4) 

A3 — i 

As before [H], we have carried out calculations using the fixed-nuclei approximation 
[22| [23] , taking the internuclear separation to be at the equilibrium value, i? = 1.4 au. 
We have labelled the electrons as particles 1 and 2, taking the positron to be particle 3. 
The position vector, r^, of each lepton is described by the prolate spheroidal coordinates 
[23] {Xj, j G {1,2,3}. These coordinates are defined implicitly in terms of the 



Cartesian coordinates, {xj,yj,Zj), as 

= [(\' - 1) cos (0,), (5) 

= [(\' - 1) sin (0,), (6) 

= l^^jf^j- (7) 



The functions S and C represent, respectively, the incident and scattered positrons 
asymptotically far from the target. The shielding parameter, 7, ensures the regularity 
of C at the origin and is taken to have the value 7 = 0.75. The constant, c, is defined to 
be c = kR/2, k being the magnitude of the positron momentum in atomic units. is 
a normalization constant and can here be regarded as arbitrary. The unknowns, at and 
{po, . . . ,Pm}, are constants to be determined. The inclusion of the parameter, r, in \l/t 
is a generalization of the Kohn method due to Kato [251 ES]. This parameter is of only 
minor physical significance, playing the role of an additive phase factor in the part of 
the wavefunction representing the incident and scattered positrons asymptotically far 
from the target. However, at each value of k, the value of r can be varied to avoid 
spurious singularities in the Kohn calculations. Away from the spurious singularities, 
for an accurate trial wavefunction we can expect the variation in the calculated values 
of r/v over r to be small. In the original application of the Kohn method [1], only 
wavefunctions corresponding to r = were considered. 

The function, ipQ, is an approximation to the ground state wavefunction of the 
unperturbed hydrogen molecule and is determined by the Rayleigh-Ritz variational 
method [27] . In the calculations presented here, we have taken ipQ to be the target 
wavefunction described in detail in another of our previous calculations [10] , accounting 
for 96.8% of the correlation energy [27] of H2. The function, 

A^ 

XO = T 7C0S[C(A3 - 1)] 

A3 — i 

X {1 - exp [-7 (A3 - 1)]} exp [-7 (A3 - 1)] , (8) 

is the same as has been used in our earlier calculations [9], |T0] and was introduced 
first by Massey and Ridley [3]. The remaining short-range correlation functions. 
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Q = {xi;---)Xm}; allow for the description of direct electron-positron and electron- 
electron interactions. Here, we have used the same set of M = 279 correlation functions 
described in detail in equations (5-8) of [lO]. The general form of each function, Xj, is 

Xi = /i(ri,r2,r3)exp[-/5(Ai + A2)-aA3] (1 < « < M) , (9) 

where each /j (ri,r2,r3) is symmetric in the coordinates of the electrons. They are a 
mixture of separable correlation functions and Hylleraas-type functions [28] containing 
the electron-positron distance as a linear factor. As discussed previously [10], the 
Hylleraas-type functions in particular allow for high accuracy of results away from 
anomalous singularities. Unless otherwise noted, we have here chosen values of a = 0.6 
and /5 = 1.0, rather than the values of a = 0.3 and (3 = 0.7 used earlier [10]. This choice 
of nonlinear parameters highlights the interesting aspects of Schwartz-type anomalies 
more clearly. 

In our application of the Kohn variational principle [HI El], the functional 

J [vft] = tan (r],-T + c) = a,- -^J^(v[/,, v]/,) (10) 

is made stationary with respect to variations in a^ and {po, . . . ,Pm}- Here, we have 
denoted by (\E't, ^t) the integral, (\&t| (-^ ~ 1^*)' where H is the nonrelativistic 
Hamiltonian for the scattering system and E is the sum of the positron kinetic energy 
and the ground state energy expectation value of ipQ. The integral is evaluated over the 
configuration space of the positron and the two electrons. We will, henceforth, use this 
notation more generally to denote by (X, Y) integrals of the form {X\ (^H — \Y). 
The stationary principle imposed upon (fTOl) leads to the linear system of equations 

Ax = -b, (11) 

where 
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Solving (ITT!) determines the values of at and {pq, . . . ,Pm}, allowing {"^t, and, hence, 
?7v to be calculated via (|T0|) . However, as has been discussed extensively (see, for 
example, P [29l [30], [3l]), the particular form of the functions, fi (ri, r2, rs), used in our 
calculations does not, in general, permit analytic evaluation of the integrals comprising 
the matrix elements of A and b. Sophisticated methods to determine these integrals 
numerically have been developed [29l [30} [31]. However, the numerical approaches can 
give only accurate approximations to the exact values of the integrals, so that small 
errors in determining the elements of A and b are unavoidable. 

2.2. Singularities 

Singularities in our generalized Kohn calculations arise from zeros of det (A), the 
determinant of f[T2|) . Under these circumstances, the linear system f[TT|) has no unique 
solution. Close to these singularities, it is well known [TTl [T2| [T3| [TH [15] that values of 
?7v obtained by solving ( ITT]) can be anomalous; small errors in the elements of A or 6 can 
correspond to large errors in the solution, x, particularly when A is close to singularity 
in a sense that can be defined formally in terms of the condition number [321 [33] of A. 
A more detailed discussion of the condition number will be given in section 13.21 

It is appropriate at this point to define a convention that we will adopt in our 
discussion of singularities in the generalized Kohn method. The type of spurious 
singularities mentioned by Schwartz [TTl [T2] here correspond to zeros in det (A) for 
the particular case when r = 0. We will, however, find it convenient to label as 
Schwartz singularities those zeros of det (A) occurring at any G [0, vr) which give rise to 
anomalous behaviour in the calculations of rj^ (t) when r is near Tg. This is an important 
clarification for the following reason: we claim that, because of our inclusion of r in \E't; 
there exist zeros of det (A) which are not spurious and which do not correspond to 
anomalous behaviour in the values of rj^. We will refer to such singularities as anomaly- 
free singularities. 

To understand how anomaly-free singularities might arise, it is helpful to consider 
the component, \l/o, of the exact scattering wave function, corresponding to the lowest 
partial wave. \E'o can be expanded as 



where ifj is the exact ground state target wavefunction and the complete set of correlation 
functions, {Ci}, describes exactly the leptonic interactions at short-range. As noted by 
Takatsuka and Fueno [I6] in their Kohn calculations for single channel scattering, the 
exact phase shift, rjQ, determined by is independent of the choice of r in f[T5l) . As a 
result, there is precisely one value, tq G [0, vr), at each positron energy such that 



oo 




(15) 



i=l 



rio-To + c = ±n- 



(16) 
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for some odd value of n > 0, where either +n or —n is chosen to keep ?7o G (— 7r/2, vr/2]. 
The value of cot {rjo — t + c) will then pass continuously through zero as r passes through 
To- 

Returning to the trial wavefunction ([1]), for nonsingular A it can be shown using 
(fTOl) and (HI]) that 

, TiN^R^k f det (A) \ 

- ^ + ^) = U-det(A)(5^o,^^c)j ' ^''^ 

where we have defined 

X = [adj (A)] &, (18) 

adj (A) being the adjugate matrix of (fT2i) . We note that adj (A) exists even if A 
is singular. In section 13.11 we will provide numerical evidence that the limit of 
cot (?7v — r + c) as det {A) exists and is equal to zero. We can, therefore, assert a 
correspondence between zeros of cot {jj^ — t + c) and zeros of det {A) . 

Suppose that, at each k, there are values of r making A singular, which we 
will denote by = {ri'^'^'', ri'^'^^ . . . , ri'^'"^'''*}. If our previous assertion is correct, each 
element of will correspond to a zero of cot (77^ — r + c). We can reasonably claim that 
if \E't is, in some sense, sufficiently accurate, precisely one of these zeros will correspond 
to the zero of cot {rjo — t + c) necessarily found at Tq for ^'q- Assuming that this is 
the case, we will denote by fs the element of corresponding to tq. Values of rj^ 
determined in the generalized Kohn method should then vary slowly and smoothly with 
r as it passes through fg. Indeed, as r ^ fg from either side, we would expect the values 
of ?7v calculated by solving ( |TT1) to converge to the value, r/v, determined directly from 

V^,-fs + c= ±n-, (19) 

where n is again chosen so that fj^ G (— 7r/2,7r/2]. 

In the following section we will present results of generalized Kohn calculations 
exhibiting anomalous behaviour due to Schwartz singularities and, further, demonstrate 
empirically that the anomaly-free singularities do exist and that values of fg can be 
found. At each k, choosing r = % (k) then defines an optimization of r that will be seen 
to avoid anomalies in rj^ due to Schwartz singularities. 



3. Results 



3.1. Calculations of phase shift 

In our generalized Kohn calculations, we have obtained values of rj^ G (— 7r/2,7r/2] 
using ffTOj) and f|TT]) . for a range of positron momenta. Spurious singularities have been 
accounted for by performing calculations over p different values of r equidistant in the 
range r G [0, vr). For the results presented here, we have taken p = 1001. Calculations for 
a large number of r values can be carried out with minimal additional computational 
effort, as it can be shown that the matrix elements of A and b for any r are readily 
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Figure 1. Values of rj^ (r) at fc = 0.2. 




2.868 2.869 2.87 2.871 2.872 



T 

Figure 2. The behaviour of cot (t^v — t + c) at fc = 0.2 for values of t either side of a 
singularity. 



available from the elements of A (r = 0) and 6 (r = 0) via an orthogonal transformation. 
It is helpful to carry out calculations for large p as it allows a detailed examination of 
the behaviour of rj^ very close to Schwartz singularities to be made. Values of ?7v over 
the p values of r are given in figure [1] for k = 0.2, corresponding to a positron energy 
of 0.54 eV. Anomalous results due to a Schwartz singularity are clearly evident around 
T ~ 2.87. We have indicated the value of r giving rise to the singularity by a dashed 
line. Away from this value of r, the variation in rj^ is small. 

In figure [2] we have studied the behaviour of cot (?7v — t + c) close to the singularity 
at r ~ 2.87, by calculating cot (77^ — r + c) for 101 values of r equidistant in the range 
T G [2.868, 2.872]. We have again indicated the position of the singularity in this figure 
by a dashed line. The results shown in the figure suggest that cot (77^ — r + c) converges 
smoothly to zero as r — > Tg from either side, supporting the assertion made in section 
12.21 regarding the correspondence of the zeros in cot {f]^ — t + c) and det (A) . We have 
found that behaviour of the type shown in figure [2] is a general feature of the calculation. 
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X 



Figure 3. Values of dot {A) at fc = 0.2 for < r < tt. 

Using ([2]) and f|T2l) . it is straightforward to show that 

det (A) = A (k) sin^ (r) + B (k) sin (r) cos (r) + C (k) cos^ (r) , (20) 

where A{k), B{k) and C (/c) are constants with respect to variations in r. For a given 
positron momentum, the constants, A, B and C, can be determined by calculating 
det {A) directly from ( fT2|) at particular values of r. Strictly speaking, in our calculations 
we have evaluated det (^A^ , A being the approximation to A whose elements have been 
determined using numerical integration. We will assume that the values of A, B and 
C are not unduly sensitive to small changes in the elements of A and henceforth take 
det {A) and det {^A^ to be essentially equivalent. 

At each /c, provided that ^ 7^ 0, the values, Tg, making A singular can be found by 
solving the quadratic equation in tan (ts), 

^tan2(rs) + i3tan(rs) +C = 0. (21) 

If only r is varied, and unless A = B = C = there will be no more than two zeros of 
det [A) in the range r G [0, vr). Figure [3] shows det {A) as a function of r at /c = 0.2. 
The scale on the vertical axis is unimportant, since the value of det {A) at each k can 
be made arbitrarily large or small by a choice of the normalization constant, A^. The 
result of interest in the figure is that it indicates two values of r at which A is singular. 
The anomalous behaviour in figure [1] corresponds directly to the singularity observed at 
r ~ 2.87 in figure O However, there are no anomalies in figure [U corresponding to the 
singularity at r ~ 1.90 in figure [3l suggesting that this singularity is of the anomaly- free 
type described in the previous section. 

We have examined this phenomenon at other values of the positron momentum. 
Figure H] indicates the roots of fl2T|) for 100 different positron momenta equidistant in 
the range 0.01 < k < 1, corresponding to a positron energy range from 1.36 meV to 
13.6 eV. For the majority of positron momenta considered here. A, B and C are such 
that there are two values of Ts at each fc, the exceptions being k = 0.65 and k = 0.66, for 
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12 3 
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S 

Figure 4. The zeros of det (A) for 0.01 < fc < 1. 

which we have found no real- valued solutions of ( 12T]) . It is apparent from figure H] that 
the roots of fl2Tl) lie in two families of curves. The first family spans the entire range, 
Ts G [0, it), for 0.01 < k < 1. The second family is confined to values of Tg in the range 
Ts G [1.5,2.1] for all positron momenta considered here. For almost every k where real 
roots exist, there is precisely one root corresponding to each family. The exception is 
the result ai k = 0.71, where there is an irregularity in the otherwise smooth behaviour 
of Ts over k. We will discuss this phenomenon in more detail in section 13. 2[ 

To illustrate anomalies in rjy at different values of fc, it is convenient to define the 
function, 

A{k,T) = \vAk,r)-{Vy){k)\, (22) 

where, at each k, (rj^) is the median value of rj^ (r) evaluated over the p values of r. 
A {k, t) measures the degree to which a given ri^ {k, r) can be considered anomalous. 
The values of A [k, r) are shown in figure [5l For clarity, we have included results only 
for 50 values of k equidistant in the range 0.02 < A; < 1, rather than the 100 values 
used for figure HI The omission of the results for k = 0.71 in this figure also allows us to 
delay until section 13.21 the discussion of the atypical singularity observed at this value 
of k in figure HI 

It is clear from figure [5] that anomalies are observed corresponding to only the first 
of the two families of curves identified from figure HI This is strong evidence that the 
curve for which no anomalies are observed comprises legitimately occurring singularities. 
It is interesting to note that the size of the anomalies due to the Schwartz singularities 
become noticeably smaller in figure [5] as they tend to coincide with the apparently 
anomaly-free singularities. Denoting by ri^'' (k) the values of r describing the anomaly- 
free curve, at each k we expect that, as r — > Ts^\ the values of rj^ calculated over r will 
converge smoothly to the value, r)v, determined directly from (flQl) and taking fg = ri^\ 
We have found in our calculations that this is indeed the case, and an example is shown 
in figure [6] for k = 0.2. 
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Figure 5. Values of A {k, r) for 0.02 < fc < 1. 
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Figure 6. Convergence of [ — ]77v (t) to [+]77v as r — > fg at fc 0.2. 

We can reasonably conclude that we have encountered anomaly-free singularities. 
Before we can develop an optimization for r based on these singularities, however, there 
are two outstanding issues to be addressed. Firstly, we have already noted that no 
real- valued solutions of fl2ip were found at k = 0.65 and k = 0.66. This is potentially 
problematic as, for a sufficiently accurate trial function, we should expect at least one 
real root of ( l2Ti) at each k, corresponding to an anomaly-free singularity. However, 
inspection of figure H] shows that the solutions of (!2T]) are close together in the regions 
either side of /c = 0.65 and k = 0.66. Were the two roots to coincide at some k, then 
= 4AC. Near a point of coincidence, ~ 4AC and small errors in the values of 
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Figure 7. Values of 3 [ts] for 0.64 < fc < 0.67. 



A, B and C could erroneously give rise to — AAC) < 0. The values of r solving 
fl2T]) at A; = 0.65 and k = 0.66 were found to be, respectively, Ts ~ 2.02 ± O.Oli and 
Ts ~ 1.98 ± 0.03i. In both cases, the fact that ^3 [ts] <^ 3? [ts] suggests that singularities 
do genuinely exist for r G M at these values of k, but small errors in our calculations of 
A, B and C due to inexact numerical integration have prevented us from finding them. 
Having investigated this problem in more detail, in figure [7] we show the calculated 
values of [ts] for 31 values of k equidistant in the range 0.64 < k < 0.67. There is a 
clearly defined region of k where no real roots of ( 12T|) have been found. The smoothness 
of 53 [ts] over k in this region does not necessarily preclude the notion that the failure 
to find real-valued solutions is due to small numerical errors in our calculations. It is 
conceivable that inaccuracies in the calculated values of A, B and C could also arise from 
systematic errors in the algorithm [31] used to calculate the determinants. Nevertheless, 
the results illustrated in figure [7] are interesting; their exact origin may be speculated 
upon and will remain a subject of our ongoing investigations. 

The second difficulty concerns the choice of fs from the two available solutions of 
f l2T|) . A method is needed for identifying at each k the root of fl2Tl) corresponding to a 
legitimate singularity. This can easily be achieved by inspecting values of r/v at values 
of r either side of each singularity, although this approach is not ideal as it requires 
solutions of (ITT]) to be found. In practice, at each k it should be possible to determine 
by inspection which of the two phase shifts is anomaly-free by examining corresponding 
results at singularities for nearby values of k. For example, figure H] clearly shows that 
only one curve in the (r, k) plane corresponds to a physically acceptable variation of 
phase shift over k. 

With these considerations in mind, we claim that choosing r = fg at each k defines a 
consistent optimization that can be used to avoid anomalies due to Schwartz singularities 
appearing at other values of r. To evaluate the success of this approach, we have found 
it helpful to consider an alternative optimization of r. For the calculations of rj^ carried 
out with p = 1001, choosing the value of r at each k giving rise to the median phase 
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Figure 8. A comparison of optimization schemes for r, [ ] (77v) (fc) and [x] fi^ (fc). 

shift, (?7v), should also mitigate anomalous behaviour. In figure [H] we have compared 
results for (?7v) (k) and fj^ (k) for momenta in the range 0.01 < A; < 1. Here, of the two 
candidates for fjy corresponding to the two singularities, at each k we have chosen the 
one whose absolute value is closest to |(77v)|- For clarity, we have included in the figure 
values of fj^ (k) for only 50 values of k equidistant in the range 0.01 < k < 0.99. 

Both optimization schemes successfully avoid anomalous behaviour at most values 
of k and there is good agreement between the two sets of results at all momenta. 
However, the intriguing feature of the figure is the anomaly appearing in both sets 
of results at /c ~ 0.71. For (77^) (k), we have shown this anomaly in greater detail by 
including in the figure results of a further 100 Kohn calculations for momenta equidistant 
in the range 0.7 < A; < 0.72, although it is practical to consider henceforth only the 
anomalous behaviour occurring precisely a.t k = 0.71. We believe the anomaly shown in 
figure [8] is of a different type to those shown in figure [5l which are due to singularities 
found by varying only r at a given k. In the following section we will examine the 
circumstances under which persistent anomalies of the kind shown in figure M could 
appear, before going on to discuss methods designed to avoid them. 

3.2. Persistent anomalous behaviour 

It is often claimed that anomalous results observed in the region of singularities arise 
from A having a determinant close to zero. Statements of this kind can be misleading, 
as the determinant of any nonsingular A can be made arbitrarily close to zero by an 
appropriate scalar multiplication, without altering the sensitivity of the solution, x, to 
small errors in the elements of A ot b. A better measure for identifying regions where 
anomalies may occur is the condition number, k{A), defined for nonsingular A. The 
condition number is independent of the normalization constant, A^, and for the linear 
system (fTTl) is defined as the maximum ratio of the relative error in x and the relative 
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error in b. Formally, it can be shown that 

k{A) =\\ A\\\\ A~'^ \\, (23) 

with respect to some matrix norm, || A \\ [32]. The value of k, (A) is dependent upon the 
choice of norm. In our calculations, we have considered the matrix 1-norm [32] of A, 



M+2 

II A II 1= max > loiJ, (24) 

l<i<(A/+2) ^ ' ^' 

where is the element in the i^^ row and j*^ column of A. In what follows, the 
particular choice of the 1-norm in our calculations will implicitly be assumed. 

A matrix with a large condition number is said to be ill-conditioned, and the solution 
of the corresponding linear system may not be reliable if the elements of A and b are 
not known exactly. For any invertible A, the condition number can be used to formalize 
the definition of closeness to singularity in the following way. If AA is defined to be any 
matrix such that A-l- is singular, then the relative distance to singularity, A (A), for 
A, is defined [32] to be 

A (A) = min (^^^ ■ det {A + AA) = 0^ . (25) 

This definition holds for any consistent norm. Further, if A (A) and k (A) are evaluated 
using the same choice of norm, it can be shown [32] that 

A{A) = [k{A)]-\ (26) 

In section 13.11 we noted that (!2T!) has, in general, no more than two zeros if only 
variations of r G [0, vr) are considered. However, ifA = B = C = then det (A) 
is identically zero independently of r and no consistent value of phase shift can be 
calculated, either by solving (fTTl) or directly from ( |T9l) . There is no obvious physical 
reason why this circumstance should arise at any /c G M. However, it is conceivable that 
A, B and C could coincident ally be close to zero, in some sense, over a narrow range of 
k. Small errors in the evaluation of A, B and C could then give rise to both unreliable 
solutions of fl2T|) and persistent anomalies in the calculation of rj^ due to ill-conditioning 
in the Kohn equations ( |TT|) . To see how the latter case arises, using ( l20l) we note that 



C 
A 

B 



det [A (r 
A 



det 
2 det 



0)], 

TT 



A-C, 



(27) 
(28) 

(29) 



so that the notion of A, B and C being close to zero is immediately formalized in terms 
of A (A) at r = 0, r = 7r/4 and r = tt/2. 

In our calculations, we have used a numerical algorithm [35] to calculate Aa(A), 
an estimate of A (A). Values of Aa (A) at r = 0, r = 7r/4 and t = ir/2 are shown in 
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k 

Figure 9. Values of Aa (A) over fc for [ ] r = 0, [ — • — ] t ~ j and [ ] t ~ ^. 

figure [9] for 100 values of k in the range 0.7 < k < 0.72. Values of Aa (A) for r = 
and r = 7r/4 are anomalously small at A; ~ 0.71, with the values for (A) at r = 7r/2 
also passing through a clear minimum at k ~ 0.711. Ordinarily, we would not expect 
ill-conditioning to occur over a very broad range of r at a given k. At k = 0.71, the 
small values of Aa (A) atr = 0,r = 7r/4 and t = 71/2 therefore point to a manifestation 
of ill-conditioning which is unusually widespread in r. In fact, in our calculations we 
have failed to find any value of r G [0, vr) such that A is sufficiently well-conditioned to 
avoid anomalous results at = 0.71. We have also confirmed that A, B and C all pass 
through zero at least once between k = 0.71 and k = 0.7104. 

We can conclude that both optimization schemes developed in section 13.11 
successfully avoid anomalies due to Schwartz singularities whose existence depends only 
upon the choice of r. However, at certain values of k, A can become close to singularities 
whose existence is independent of r, resulting in anomalies in the calculation of t]^ that 
persist even after r has been optimized. In attempting to address this problem, a 
number of methods are available. In principle, anomalous behaviour can be reduced 
dramatically by allowing the Kohn trial wavefunction to be complex- valued [U [17] . An 
alternative approach for avoiding anomalies persistent in r is to make a small change in 
some other parameter of the trial wavefunction. We will explore both of these techniques 
in the following section. 

3. 3. The complex Kohn method 

The complex Kohn method is an extension of the original variational approach in which 
the boundary conditions of the trial wavefunction are complex. It was originally believed 
[H [17] that this method was anomaly-free, although anomalies were subsequently 
reported by Lucchese [18]. For our complex Kohn calculations on (e"*" — H2) scattering. 



o 

6 



< 



2 
n 
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we have used a trial wavef unction, of the form 

AI 



^[={S + a[f + p'oxo) + 5^ p\Xii (30) 



=1 



where 



T = 5 + iC, (31) 



the functions ipQ and VL = {xi, ■ ■ ■ ,Xm} being the same as in ([T]). The unknowns a[ 
and {pq, . . . ,p\f} will not, in general, be real. Application of the variational principle 
to (130|) leads to a matrix equation analogous to (fTTj) . 



A'x' = -b' (32) 

where A' and 6' are identical to A and 6, but for the function, T, replacing C in f|T2l) 
and ( fT3l) . The determinant, det (A')' conveniently reduces to 



where 



det {A') = V (k) e-^'\ (33) 



V={A-C)-iB (34) 



is a complex constant with respect to variations in r. The values of det [A') then 
describes a circle of radius |D| in the complex plane for variations of r G [0, vr). Hence, 
singularities are obtained only if both the real and imaginary parts of V are zero; they 
can neither be located nor avoided by varying only r. We would therefore expect 
anomalous results due to singularities arising from the choice of r at a given k to be 
eliminated in the complex Kohn method. However, from flM|) we would also expect T) to 
be close to zero at /c = 0.71, in the same sense that A, B and C have already been seen 
to be close to zero at this value oi k. It is therefore likely that the anomalies already 
seen to occur due to relationships between A., B and C rather than the choice of r, will 
persist even in the complex Kohn method. 

We have obtained values of rj^ using the trial function, We have found that the 
differences in the calculated values of rj^ at different values of r G [0, tt) are negligible, 
for all positron momenta considered here. Without loss of generality, we can regard the 
complex Kohn calculation as effectively independent of r and choose r = for simplicity. 
In figure fTOl we have compared results for rj^ {k, r = 0) obtained with the trial function, 
'^[, with the results for (r/v) (k) obtained in section \3A] with 

The results of the two methods are essentially equivalent at all positron momenta, 
with the optimized results for \Eft differing from the results for \E^( by no more than 0.1% 
at each k. As expected, the complex Kohn method automatically avoids anomalies at 
most values of k without the need for an optimization of r. Nevertheless, the use of 
the complex trial function has predictably failed to remove the persistent anomaly at 
k = 0.71. In figure fTTl we verify that Aa {A') is anomalously small at /c ~ 0.71 for r = 0. 
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Figure 10. A comparison of [ ] (?7v) (fc), determined using the trial function, 

^t, with [x] ?7v (fc,T — 0), evaluated using the complex Kohn method with the trial 



function, ^![. 
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Figure 11. Values of Aa {A') over fc for r = 0. 



We have found that differences between the results shown in figure [H] and values of 
Aa {A') calculated at other values of r are negligible. 

Having failed to find a systematic remedy for the persistent anomalous behaviour, 
we consider a more ad hoc approach. It should be possible to avoid any Schwartz-type 
anomaly by some variation of parameters in the trial wavefunction. We have found that 
variations in r are not always successful, but other candidates exist. In our complex 
Kohn calculations, we have varied the values of a and in ([9]), fixing r = 0. Recall that 
the values of the these parameters have so far remained fixed at a = 0.6 and (3 = 1.0. 
We now consider the results of Kohn calculations carried out for 31 different values of 
a in the range a G [0.59, 0.605] and 61 different values of jS in the range /5 G [0.65, 1.25]. 

To illustrate persistent anomalous behaviour, it is helpful to define a function 
analogous to fl22|) . 



A' {a, (3) = |?7v (a, P) - (") 



(35) 
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Figure 12. Values of A' {a, /?) at fc = 0.71. 

0.74 I ■ ■ — 




0.68' ■ ■ ■ ■ 1 
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y 

Figure 13. Values of rj^ (7) at fc = 0.71. 

where, for each of the values of a considered, fj^ (a) is the median value of r/v evaluated 
across the range of values of /3. Values of A' («, f3) are shown in figure [T2| from which it 
is clear that persistent anomalies appear distributed about a curve in the {a, (3) plane. 
For values of a and (3 away from this curve, the calculations are free of anomalies. 
Hence, a small change in the values of a or /5 can indeed be shown to successfully avoid 
persistent anomalous behaviour. 

Finally, we consider briefly that the shielding parameter, 7, in (jlj) and ([8]) might 
also be varied in an effort to avoid anomalous behaviour. Values of t]^ a.t k = 0.71, 
a = 0.6 and /3 = 1.0 for 0.5 < 7 < 1.0 are shown in figure [131 It is apparent that 
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small changes in the value of 7 have relatively little effect on the persistent anomaly at 
k = 0.71. This is not unexpected, being consistent with the findings of Lucchese |18j . 
who investigated the effect of varying a parameter analogous to 7 in his model potential 
calculations. He noted that singularities due to the choice of 7 occurred when the values 
of 7 and a were not similar, most typically when 7 <C a. With this in mind, and from 
inspection of figures [12] and [13], we can conclude that the anomaly observed in figure [10] 
is due primarily to the choices of a and (3 rather than the choice of 7. 

4. Concluding remarks 

We have carried out a thorough examination of singularities and related anomalous 
behaviour in generalized Kohn calculations for {e~^ — H2) scattering. We have argued 
that singularities do not always occur spuriously and that variational calculations of the 
scattering phase shift can be anomaly-free at these singularities. Subsequently, we have 
developed an optimization scheme for choosing a free parameter of the trial wavefunction 
allowing anomaly-free values of the phase shift to be determined without the need to 
solve the linear system of equations derived from the Kohn variational principle. This 
approach has been seen to be largely successful, giving phase shifts in close agreement 
with those determined by a conventional generalization of the Kohn method, as well as 
those obtained with the complex Kohn method. 

Persistent anomalies in both sets of calculations have been identified and attributed 
to singularities that cannot be avoided with any choice of the parameter, r. Further, 
we have found that our implementation of the complex Kohn method is susceptible to 
the same behaviour. We have demonstrated, however, that persistent anomalies can 
be avoided by small changes in the nonlinear parameters of the short-range correlation 
functions. Hence, by studying the behaviour of A, B and C over we can predict 
the appearance of persistent anomalous behaviour quantitatively and avoid it by an 
appropriate change in a or (3. 
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